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Abstract. Cosmic voids in the large-scale structure of the Universe affect the peculiar mo¬ 
tions of objects in their vicinity. Although these motions are difficult to observe directly, the 
clustering pattern of their surrounding tracers in redshift space is influenced in a unique way. 
This allows to investigate the interplay between densities and velocities around voids, which 
is solely dictated by the laws of gravity. With the help of V-body simulations and derived 
mock-galaxy catalogs we calculate the average density fluctuations around voids identified 
with a watershed algorithm in redshift space and compare the results with the expectation 
from general relativity and the ACDM model. We find linear theory to work remarkably well 
in describing the dynamics of voids. Adopting a Bayesian inference framework, we explore the 
full posterior of our model parameters and forecast the achievable accuracy on measurements 
of the growth rate of structure and the geometric distortion through the Alcock-Paczyhski ef¬ 
fect. Systematic errors in the latter are reduced from ^ 15% to ^ 5% when peculiar velocities 
are taken into account. The relative parameter uncertainties in galaxy surveys with number 
densities comparable to the SDSS MAIN (CMASS) sample probing a volume of l/i“^Gpc^ 
yield cr//6/(//5) ^ 2% (20%) and ctdaH/DaH ^ 0.2% (2%), respectively. At this level of 
precision the linear-theory model becomes systematics dominated, with parameter biases that 
fall beyond these values. Nevertheless, the presented method is highly model independent; 
its viability lies in the underlying assumption of statistical isotropy of the Universe. 
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1 Introduction 

With the advent of modern galaxy redshift surveys that map out significant fractions of the 
observable Universe in ever more detail (e.g., [1-7]), its least luminous and yet most extended 
constituents have drawn an increasing amount of attention recently: cosmic voids. In fact, 
cosmic voids are not only interesting objects in their own right. As they take up most of 
space in the Universe they play a major role in the formation of the web-like pattern of 
its large-scale structure. Despite its complexity, this large-scale structure contains a wealth 
of information on the fundamental properties of the Universe, such as its initial conditions, 
its matter and energy content, or the geometry of spacetime on cosmological scales, for 
example. Yet not only cosmology, but the very nature of gravity itself can be investigated 
with large-scale structures, as it is gravitation that gives rise to structure formation in the 
first place. This happens through gravitational collapse of initially over-dense regions in the 
mass distribution into filaments, sheets and clusters. Thereby, the constituents of this mass 
distribution acquire increasing velocities over time, yielding a complex dynamical behavior 
in 6-dimensional phase space. However, the dynamics of voids remains rather simple, as they 
are mainly composed of single-stream flows [8-10]. 

Unfortunately, a dominant fraction of the mass distribution in the Universe is made 
up of invisible dark matter and we can only observe the tracers of this mass distribution. 
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such as galaxies or Hydrogen gas clouds,^ for example. The statistical spatial distribution 
of these tracers is typically different from the one of the dark matter, which goes under the 
name of bias. Therefore, in order to reconstruct the mass distribution of large-scale structure 
one has to understand the relation between tracers and the dark matter, for instance with 
the help of theoretical models or numerical simulations. Furthermore, this reconstruction is 
complicated by peculiar motions. The peculiar velocity of a tracer is difficult to determine 
observationally, and is typically not known for distant objects. However, it affects the redshift 
of each tracer via the Doppler effect in addition to the redshift caused by the Hubble ffow. 
Since the total redshift is used to determine each tracer’s line-of-sight distance and the two 
sources of redshift cannot be distinguished, the apparent distribution of tracers in redshift 
space (as opposed to real space) is skewed, an effect referred to as redshift-space distortion. 
On large scales this results in elongated structures perpendicular to the line of sight due to 
the coherent streaming motions of tracers towards overdensities in the matter distribution, 
the so-called Kaiser effect [12]. In contrast, on small scales structures preferentially stretch 
along the line of sight, which is commonly referred to as the Finger-of-God (FoG) effect. It 
is caused by nonlinear virial motions of tracers within gravitationally collapsed objects, such 
as dark matter halos. 

Similarly, incorrect assumptions about the expansion history and the geometry of the 
Universe can result in a skewed distribution of tracers. The conversion between redshifts and 
angles on the sky into Cartesian coordinates requires knowledge of the spacetime metric that 
describes the Universe, which in turn depends on the latter’s matter and energy content. Any 
deviation from the true underlying cosmological parameters results in a distortion of tracer 
coordinates, which is known as the Alcock-Paczyhski (AP) effect [13]. In cases where the 
underlying statistical arrangement of tracers is known, it can be used to constrain cosmology. 

Luckily, a tracer’s angular position on the sky is not affected by redshift-space distor¬ 
tions. Therefore, by statistically comparing its line-of-sight and angular clustering patterns, 
one can disentangle the impact of peculiar motions and geometric distortions on the coor¬ 
dinate distribution of a tracer. This method is particularly powerful for cosmic voids, as 
on average they exhibit spherical symmetry in real space, which allows to assign any type 
of anisotropy to a skewed coordinate system in redshift space [14-18]. Vice versa, these 
anisotropies can be used to infer the peculiar velocities of tracers around voids [19-22], 
opening up the possibility to determine the strength of gravity on void scales. The main 
underlying assumption to this method is the cosmological principle, i.e. statistical homo¬ 
geneity and isotropy of the Universe. Observational evidence for the anisotropic nature of 
the void-galaxy cross-correlation function has been provided by Paz et al. [21] using SDSS 
data, and Micheletti et al. [23] in the VIPERS survey. 

While the investigation of redshift-space distortions between pairs of galaxies has re¬ 
ceived enormous attention in the literature, both theoretical and observational (see e.g., [24- 
40] and [41-55], respectively), comparable studies concerning voids remain rare [18, 19, 21, 
23] . The aim of this paper is to assess the feasibility of measuring redshift-space distortions 
around voids using galaxies from a realistic mock galaxy catalog, and to compare this data 
with a physically motivated model. The values for the parameters that go into this model 
will be investigated in detail, including their uncertainty, in order to forecast the potential 
to test the predictions of General Relativity and the current standard model of cosmology 
with the help of modern surveys. 

^See ref. [11] for an interesting proposal to identify voids in the Lyman-alpha forest. 
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2 Theoretical prerequisites 


2.1 Tracer coordinates in redshift space 

Let us denote the coordinates of a tracer in real space by the vector x. We can only observe 
its angular position on the sky and its redshift 2 ;. If this redshift were only caused by Hubble 
expansion, we could compute the comoving line-of-sight distance to the tracer as 


where c is the speed of light and H{z) describes the Hubble expansion as a function of 
redshift. The tracer’s peculiar velocity component along the line of sight causes a small 
additional Doppler shift of magnitude which is indistinguishable from the cosmological 
redshift. Thus, the apparent comoving line-of-sight coordinate of a tracer in redshift space 
becomes 


X\\ = 


I 


z+^(l+z) 


dz' ~ xii + (1 + z) . 


( 2 . 2 ) 


/o H{z') " H{z) 

Furthermore, in cosmological redshift surveys the angles d between tracers observed on the 
sky are converted to comoving distances x_\_ using the angular diameter distance Da{z), 


x± = . 


(2.3) 


In a Friedmann-Lemaitre-Robertson-Walker metric of curvature k the comoving angular di¬ 
ameter distance is given by 

It is sensitive to the curvature parameter = 1 — Dm — Da, while the Hubble rate itself 
also depends on the matter and energy content in the Universe, 

H{z) = • (2-5) 

It is hence necessary to assume fiducial values for today’s Hubble constant iLo, matter fraction 
Dm, curvature parameter D^i^, and cosmological constant Da in order to construct a three- 
dimensional map of the large-scale distribution of tracers^. If these fiducial parameters do not 
coincide with the true cosmological values, or the assumed ACDM-model is incomplete, the 
distances in eqs. (2.1) and (2.3) will be incorrect. This effect is commonly referred to as the 
AP test [13], as it allows to determine the true cosmological parameter values in cases where 
the geometry of the observed structures is known, like it is the case for cosmic voids [14-18]. 
However, only the combination 

Axil 1 

Ax± Da{z)H{z) ’ 

can be constrained with this test, as this is the quantity that determines the degree of 
anisotropy in the spatial distribution of the tracers. If we define Da{z) and H{z) to adopt 
the fiducial (and possibly incorrect) values for our cosmology, 6: quantifies the deviation from 
isotropy resulting from analyzing the data with an incorrect fiducial cosmology. 

^Here we consider a cosmological constant A, but this can straightforwardly be extended to an arbitrary 
equation of state for dark energy. 
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2.2 Void-tracer cross-correlation function 

Now let us consider the average real-space radial number-density profile of tracers in distance 
r from a void center at the origin, Pvt (^)5 which can also be referred to as a one-dimensionally 
stacked void. Labeling each of all Ny void centers with index i and coordinates X^, and 
similarly each of all Nt tracers with index j and coordinates Xj, we have 


Pvtjr) 

Pt 


1 S^PvtiO _ ^ V V I r'l 

(X*-x,+r), 


(2.7) 


where the first equality describes the ensemble average over all individual void density profiles 
and the second one represents a histogram of Dirac delta functions 6^ for the tracer positions 
Xj at separation r from each void center X^. This expression can then be written as a 
convolution of the number density of void centers py with the number density of tracers pt? 

= i [ E^E^^d3x = l + ^vt(r), (2.8) 

J ^Vv iVt V J Py Pt 


where V is the total observed volume, py = Ny/V^ pt = Ni/V^ and Cvt(^) denotes the 
void-tracer cross-correlation function in real space (see also ref. [56]). 

Now, in order to obtain the corresponding redshift-space correlation function, one needs 
to consider the pairwise velocity probability distribution function 7^(v, r). This function maps 
all pairs of separation r to separation f = (ry, r±y due to their relative velocity v, such that 
the correlation function in redshift space is calculated as [56] 

1 + 1(f) = y [1 + C(0] r) • (2.9) 

According to eq. (2.2) only the magnitude of the relative velocity component along the line 
of sight V|| affects the coordinates f, which reduces eq. (2.9) to a one-dimensional integral 
via the replacements 7^(v, r) ^ and d^v d^’y. What remains to be specified is the 

explicit form of 'P('L’y,r), itself a major topic of interest in the literature (e.g., [24, 57-59]). 
Numerical studies of the pairwise velocity probability distribution function of dark matter 
particles in V-body simulations suggest an exponential form at small separations r, with a 
non-negligible skewness. Towards large separations, the skewness gets weaker and 'P('L’y,r) 
approaches a more Gaussian shape; its exponential wings however, remain [24]. 


2.3 Gaussian streaming model 

Nevertheless, for the purpose of modeling observed tracers of the dark matter, such as galax¬ 
ies, a Gaussian form for the pairwise velocity probability distribution function in many cases 
remains a reasonable approximation. Galaxies are found in dark matter halos, which move 
according to the net velocities of their constituent particles. Therefore, the small-scale non¬ 
linear dynamics of individual dark matter particles do not affect the bulk motions of halos 
and hence their hosted galaxies. This leads to a more Gaussian behavior of their velocity 
statistics, which has been confirmed with the help of simulations [26, 34, 36, 37, 60, 61]. 
Moreover, as the average tracer-density fluctuations inside and around voids are rather mod¬ 
erate [62], extreme nonlinearities (e.g. from virial motions of galaxies inside clusters) are 
less severe in this case. The fact that we model a cross-correlation between void centers and 
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their surrounding tracers also mitigates nonlinear effects, since the main velocity contribu¬ 
tion comes from the tracer, while the void center is practically at rest [63]. In contrast, in 
tracer auto-correlations the velocities at two locations are being correlated with each other, 
which effectively squares their individual nonlinearities. We will therefore assume a Gaussian 
pairwise velocity probability distribution function of the form 


V{v\ur) 




( 2 . 10 ) 


with a mean of Vy{r)^ and dispersion (T„(r). This model is referred to as the Gaussian 
streaming models and was first introduced by Fisher [64] as a description for the galaxy auto¬ 
correlation function in redshift space. Here, the comoving real-space coordinates ry and r 
can be obtained from their measured redshift-space counterparts via 


rii = rii — 


H{z, 


■(1 + ^) 


= er^ 


r = 


rn + 




( 2 . 11 ) 


Further, Vy^r) is the mean streaming velocity of tracers in distance r from the void center. 
Using the conservation of mass and following linear theory it is given by [56] 

^v(r)c.-i®®rAvn.(r), (2.12) 

O I z 

where f{z) ^ is the logarithmic growth rate of density perturbations with growth 

index 7 ~ 0.55 in the standard ACDM model [65], and 


A 


vm 



/ Pvmjq) 
V Pm 



(2.13) 


is the cumulative average mass-density fluctuation of matter around voids, compared to the 
background density p^. With the help of Wbody simulations, ref. [62] have empirically 
found a universal fitting function for the nonlinear mass-density profile of voids. 


Pvmjr) 

Pm 


- 1 = 


1 — {r/rsr 

1 +(r/rv)^ 


(2.14) 


containing four free parameters: a scale radius r^, a central under-density 5c, and two slopes 
a and /3. Here, denotes the effective void radius^ which is determined from the data itself 
as described in sec. 3.3. Given the density profile of eq. (2.14) one can integrate eq. (2.13) 
analytically, which yields 


Avm(^) = 5 c 2F1 


3 3^ fr 


- 5c ( - 


cr “h 3 


2 F 1 


(a-1-3 cr-1-3 


, (2.15) 


where 2 F 1 is the Gauss hypergeometric function. For biased tracers of the mass and following 
linear theory, one can relate eq. (2.13) to the fluctuation of tracers Avt(^) with the linear 
tracer-bias parameter 5, 


Avt(r) ~ bAy^{r) . 


( 2 . 16 ) 



















Finally, note that the velocity dispersion of eq. (2.10) depends both on magnitude and di¬ 
rection of the void-tracer separation r. It can be decomposed into its components along and 
perpendicular to the line of sight, 


(j: 


:(r) = aj{r) {r\\/rf + ^^(r) [l - (ry/r)^ 


Following ref. [66] (page 515), linear theory allows to estimate and via 

Hr) . K2{ry 
3r^ 3 J ’ 

Jsjr) _ K 2 {r) - 

2r 6r^ 3 ’ 


af(r) ^ Hlf 
4(r) ^ Hlf 


where and are defined through integrals over the correlation function, 

rr roo 

Jn{r) = dq , Kn{r) = J dq . 


(2.17) 


(2.18) 

(2.19) 


( 2 . 20 ) 


In our case we can identify ^(r) as the void-matter cross-correlation function (^vm(^) and thus, 
according to eq. (2.8), with the mass-density profile around voids, eq. (2.14). Plugged into 
eq. (2.20), this can be expressed analytically as 


Jn{r) = - 2 F 1 

n 


^ n n ^ ( r 


r y 1 
rsj a + n 


2F1 


^ a + n a + n ^ f r 


( 2 . 21 ) 


Kn{r) = (^) 


^ ( 1 


r / I /3 — n 
r \ ' 1 

r 


2F1 




(3 — a — n 


P 

2 F 1 


P 


B — a — n B — a — n 
1 - - - - 1 - 1 






ry\^' 


. ( 2 . 22 ) 


However, although linear theory is able to describe the dynamics of voids remarkably well, it 
is important to point out that it does not apply in all cases. Especially the deeply under-dense 
interior of small voids is subject to nonlinear evolution, which was shown to spoil the validity 
of eq. (2.12) at distances r < 10h~^ Mpc from the void center [62]. Moreover, the empirical 
fitting function of eq. (2.14) for the void-density profile only describes the immediate void 
environment, but not the large-scale clustering regime out to arbitrarily large scales [67, 68], 
such as the baryonic acoustic oscillation (BAO) scale, for instance. Therefore, quantities that 
depend on integrals over the correlation function out to infinity, such as Kn{r) in eq. (2.20), 
cannot be determined accurately in this case. 


3 Numerical analysis 
3.1 Simulation setup 

We analyze a large simulation that evolved N = 2048^ cold dark matter particles in a box 
of l/i“^Gpc side length with the adaptive Wbody tree-code 2HOT [69] (this simulation is 
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identical to the one used in refs. [62, 67]). It adopts the cosmological parameters of the Planck 
2013 measurement [70] with = 0.318, — 0.049, Qa = 0.682, h = 0.67, as = 0.834, 

ris = 0.962. The power spectrum of initial density perturbations is calculated with the 
Boltzmann code CLASS [71], a realization of this power spectrum is then generated with a 
modified version of 2LPTic [72] using second-order Lagrangian perturbation theory. Catalogs 
of dark matter halos and sub-halos are created using the ROCKSTAR halo finder [73]. This 
algorithm is based on an adaptive hierarchical refinement of friends-of-friends groups in phase 
space to define virialized objects. In our simulation we find about 2 x 10^ halos with masses 
ranging from ^ 2 x 10 ^^/i“^Mq to ^ 5 x 10^^/i“^Mq at redshift z = 0. 

3.2 Mock-galaxy catalogs 

In order to generate realistic mock galaxy samples we refer to a standard halo occupation 
distribution (HOD) model using the code developed by ref. [74]. It assigns central and satellite 
galaxies to each dark matter halo of mass M according to a parametrized distribution. We 
use the parametrization of refs. [75, 76], where the mean numbers of centrals and satellites 
per host halo are given by 


(7Vcen(M)) 


(A^sat(M)) 


^ 1 -h erf 

(A^cen(M)) 


/ log M — log Mn 
\ ^log M 

M-Mo^ ^ 

Ml 


(3.1) 

(3.2) 


The probability distribution of central galaxies is a nearest-integer distribution (i.e., all halos 
above a given mass threshold host a central galaxy), and satellite galaxies follow Poisson 
statistics. Central galaxies are assigned peculiar velocities of their host halo, and satellites 
are given an additional random velocity drawn from a Maxwell-Boltzmann distribution with 
a standard deviation that matches the velocity dispersion of their host halo’s dark matter 
particles. Moreover, we apply eq. (2.2) to shift all galaxies from real to redshift space. 
Without loss of generality, we pick one axis of the simulation box as the line-of-sight direction. 
We adopt the the so-called plane-parallel approximation^ where changes in the line-of-sight 
direction with the observed angles on the sky are neglected. On the scales of voids and their 
immediate surroundings, as considered in this paper, this approximation is accurate. 

We construct two different catalogs to investigate the impact of survey properties on our 
final results. For the first mock catalog we model a dense galaxy survey with the parameters 
Mmin - 2.0 X Mo 6.9 x lO^^M^Mo, Mi - 3.8 x IO^^M^Mq, ai^gM ^ 0.21, 

and cr ^ 1.12 adapted to the SDSS DR7 MAIN sample [75] at redshift z = 0. This results 
in a distribution of about 2 x 10^ galaxies of mean number density n 2.0 x 10“^/i^ Mpc“^, 
a mean separation of roughly 3.7/i“^ Mpc, and a linear bias parameter of b ^ 0.8. The 
second mock catalog is adapted to the sparser CMASS sample of the SDSS DR9 [77] with 
parameters Mmin — 1-2 x 10 ^^/i“^Mq, Mq ^ 1.2 x 10 ^^/i“^Mq, Mi 1.0 x 10 ^^/i“^Mq, 
o'logM — 0.60, and a 1.01 at redshift z = 0.5. It yields a distribution of roughly 3 x 10^ 
galaxies of mean number density h 3.0 x 10“^/i^ Mpc“^ that are separated by 15.0/i“^ Mpc 
on average, and have a linear bias parameter of b ^ 1.8. 


3.3 Void catalogs 

These mock galaxies are then utilized to generate voids with the Void IDentifieation and 
Examination toolkit VIDE [78], which in its core is based on the ZOBOV [79] algorithm. 


- 7 - 






Moreover, we use a subsampled version of our original dark matter simulation in real space 
that matches our dense mock-galaxy catalog’s sampling density of n 2.0 x Mpc“^ 

to define a dark matter void catalog as a reference. VIDE finds density minima in a Voronoi 
tessellation of the tracer particles and grows basins around them applying the watershed 
transform [80]. It gives rise to a nested hierarchy of voids and sub-voids, all of which we 
consider in our analysis. We cease to merge basins with each other if the minimum ridge 
density between them is larger than 20 % of the mean tracer density, which prevents voids 
of growing too deep into over-dense structures [79]. This value is inspired by the spherical 
evolution model [81], which predicts a minimum under-density of 0 . 2 pni for voids whose 
surrounding ridges begin to undergo shell crossing. However, this is only a reference value, as 
the watershed algorithm identifies basins irrespective of whether shell crossing has occurred. 
We further define each void center as the mean of all its member particle positions Xj, 
weighted by their individual Voronoi cell-volume Vj [14], 


. _ EjxWi 


(3.3) 


The advantage of this definition over the location of the minimum density (largest Voronoi 
cell) is its insensitivity to Poisson noise and the fact that it makes use of the watershed 
geometry which defines the void. Further, an effective void radius can be defined as the 
radius of a sphere comprising the same volume as the sum of all Voronoi cells that each void 
is composed of. 



We consider a range of Mpc ^ < 85h~^ Mpc for the voids identified in our dense 

tracer samples, and a range of 30h~^ Mpc ^ ^ 107h~^ Mpc for the ones found in our 

sparse mock catalog. This corresponds to about twice the mean particle separation as a 
lower bound, and the largest void found as an upper bound, respectively. The lower bound is 
imposed to avoid contamination of the sample by misidentified voids arising through either 
random Poisson fluctuations [79], or the effects of redshift-space distortions [82]. It amounts 
to a total number of ^ 9.9 x 10^ voids in the dark matter sample, ^ 6.4 x 10^ voids in our 
dense mock-galaxy catalog, and ^ 2.0 x 10^ in the sparse one. 


3.4 Void stacks 

As described by eq. (2.9), the tracer-density profile around voids (respectively the void-tracer 
cross-correlation function) is no longer isotropic in redshift space. The spherical symmetry of 
voids is thus reduced to a cylindrical symmetry along the line-of-sight direction. Therefore, 
it is instructive to visualize this symmetry as a two-dimensional stack with axes showing the 
void-centric distances of tracers along and perpendicular to the line of sight. To this end we 
determine the relative separation vector 


^ij = 



(3.5) 


between every void position and every tracer position Xj in units of ry, within a narrow 
range of the effective void radius with mean value fy. This is done only up to a maximum 
distance of |r^j| = 3. Further, we extract its line-of-sight component as well as its angular 




component and histogram the number of tracers in bins of width 5f|| and 5r±^ to estimate 
the void-centric tracer density according to a Poisson-counts estimator [83], 


Pvt {r\\,r_L) = 


Ei,,' 1 


4:7rf±6f±6f^^Nyr: 


(3.6) 


where Q{rij) = 0[rij — (f — 6f/2)]0[—fij -h (r + Sf/2)] combines two Heaviside step functions 
9 to define each radial bin for ry and respectively. The uncertainty in a given bin can 
then be estimated via the standard deviation in the number of tracers among all Ny voids. 
Similarly, the two components of the velocity dispersion cr||(r) and a±{r) in eq. (2.17) can be 
estimated via 




vf (fij) - v\\^±{r) 


12 


14- 




(3.7) 


except that each of them only depends on the absolute distance r to the void center and 
we weight by the Voronoi cell-volume Vj of each tracer to obtain a volumetric and hence 
unbiased estimate of these velocity statistics. The mean velocities used in eq. (3.7) are 
estimated accordingly. 




Eij ^{nj)vf^{rij)Vj 


(3.8) 


Note that we neglected the velocities of the void centers, which have been shown to be 
vanishingly small in simulations [63]. 


3.4.1 Dark matter 

Figure 1 shows two-dimensional void stacks in the dark matter for eight contiguous bins in 
effective void radius, resulting in about 12500 voids per stack. The bin edges are chosen to 
yield an equal number of voids per bin, ensuring the signal-to-noise ratio to remain similar 
across the bins. Although these dark matter stacks are not directly observable in practice, 
they shall serve as a reference to compare with our theoretical model. To this end, we 
identified physical voids in the real-space dark matter distribution, and only used the dark 
matter in redshift space for the stacking in eq. (3.6). From the figure we can distinguish two 
regimes that are affected differently by redshift-space distortions. Small and intermediate- 
size voids with Mpc exhibit a relatively high ridge that decays away towards 

decreasing angles to the line of sight (i.e. small r±). This behavior is reversed for the larger 
voids, whose ridge appears more prominent along the line of sight. These regimes may be 
attributed to the two well-known features observed in the galaxy auto-correlation function in 
redshift space: the FoG effect and the Kaiser effect [12]. The former is caused by nonlinear 
virial motions of galaxies within highly over-dense structures (such as clusters) and results in 
a stretching of contour lines on small distances along the line of sight. Because small voids 
exhibit the highest ridge densities, the FoG effect is most prominent for them and smears out 
the ridge density along the line of sight. This smearing is described by the velocity dispersion 
in eq. (2.10). 

In contrast, the Kaiser effect is caused by the coherent streaming motion of galaxies 
towards over-densities on larger scales, as described by eq. (2.12). This results in a squashing 
of contours towards the void ridges, which leads to their enhancement along the line of 
sight. Hence, the two effects act against each other and can even balance out in their 
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transition regime, which can be seen in the bottom left panel of fig. 1 with 14.9/i“^ Mpc < 
Vy < Mpc. Towards larger voids the Kaiser effect is still visible not only in the 

amplification of ridges, but also in emptying out the interior of the void, resulting in elongated 
negative contours along the line of sight (see also [19, 21, 23]). 

The influence of the AP effect on the two-dimensional void stacks is simpler to under¬ 
stand, which is why we do not show it explicitly here. According to eq. (2.6), the ratio of 
tracer coordinates along and perpendicular to the line of sight is scaled by some constant 
factor, which results in a stretching of the void stack along either ry or r±. Therefore, ge¬ 
ometric distortions due to the AP effect are expected to be far less dependent on void size 
and tracer density when compared to the dynamic distortions discussed above. 

Let us now inspect the velocity dispersion in the dark matter around these voids. The 
upper panels in fig. 4 depict its two components (j|| and a± as a function of distance r to the 
void center for the same bins in void effective radius. Both components increase rapidly up to 
the mean void radius fy of each stack, and then level off with a more moderate slope. Due to 
the rescaling of the abscissa by fy, larger voids appear to have higher velocity dispersions in 
their surrounding, but at given physical distance these values are in fact remarkably similar. 
Likewise, the differences between cr||(r) and a±{r) are quite small. In order to model these 
velocity dispersions we tried fitting eqs. (2.18) and (2.19) to our data. While this works well 
on large distances r from the void center, the velocity dispersion in the void interior is not 
accurately reproduced with these formulae, which is likely due to the mentioned discrepancy 
of linear theory on those scales and long-range correlations in the velocity statistics. Instead, 
we find that a simple Lorentzian of width oj empirically describes the observed behavior very 
well. 





1 /^ ) 

1 -h /uS^ j 


(3.9) 


Here, is simply a constant number and corresponds to the velocity dispersion at r ^ oc. 
Figure 4 depicts the best fits of eq. (3.9) to the data as solid lines. The free parameters in 
the fit are and cj, going from small to large voids they roughly range as G [230, 300] 
and uj G [1.0,0.7]fy for both cr||(r) and cr^(r). 


3.4.2 Mock galaxies 

The observationally more relevant scenario is the one where voids are identified from the 
distribution of galaxies in redshift space. Figure 2 shows the corresponding void stacks in 
our dense mock-galaxy sample at redshift z = 0 using a similar binning as before, yielding 
about 8000 voids per stack. On the one hand we notice the features in these stacks to appear 
more extreme, both the void ridges and their interiors. This can partly be attributed to 
the fact that galaxies are biased tracers of the mass, so fluctuations about the mean density 
are enhanced on void scales [84, 85]. On the other hand, this affects the redshift-space 
distortions in favor of the Kaiser effect, because galaxies are less sensitive to small-scale 
velocity dispersion in the large-scale structure as compared to dark matter particles. We 
can already observe a balance between the FoG and the Kaiser effect in the first bin at 
7Ah~^ Mpc < Ty < 8.8h~^ Mpc. For all larger voids the Kaiser effect dominates and 
causes an enhanced ridge feature along the line of sight, which can be directly related to 
the logarithmic growth rate / in eq. (2.12). For our sparser mock catalog with roughly 250 
voids per stack at redshift z = 0.5 the FoG effect is already negligible in the smallest bin, 
as apparent from fig. 3. This is partly due to the fact that we can only resolve voids with 
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Ty > 30h ^ Mpc in this sample, and partly due to the higher redshift. The Kaiser effect can 
even still be perceived for voids of effective radius Vy > lOO/i ^ Mpc. 

The velocity dispersion in the case of dense mock galaxies is shown in the middle panels 
of fig. 4. Compared to the dark matter shown in the upper panels, the behavior of cr||(r) 
and cr_\_{r) is almost identical, both in magnitude and shape of these two curves. This is 
an encouraging result, insofar as it does not suggest any velocity bias between our dense 
mock galaxies and the underlying mass distribution. Moreover, fitting eq. (3.9) to this mock- 
galaxy data yields roughly the same parameter values for ay and uu as obtained from the dark 
matter. For the velocity dispersion around voids in our sparse catalog we observe slightly 
higher amplitudes and lower widths, as can be seen in the lower panels of fig. 4. This behavior 
is expected for a more biased sample of galaxies {b 1.8), since these specifically trace out 
higher over-densities of the mass distribution that have evolved more nonlinear. 


3.5 MCMC analysis 

The question of how well our model from sec. 2.3 can reproduce the data presented in sec. 3.4 
and how well its model parameters can be constrained will be the topic of this section. For 
this purpose we implemented a Markov Chain Monte Carlo (MCMC) framework making use 
of the publicly available PyMC code package in Python [86]. Following Bayesian statistics, 
we formulate a likelihood function C that quantifies the probability of the data, given a 
model. In our case the data is given by the simulated void-tracer cross-correlation function 
Cvt(f^|| 5 f±) in redshift space, and the model is expressed by eq. (2.9) with the ingredients from 
sec. 2.3. Our model parameters can be summarized in the following parameter vector. 


e = {rs,Sc,a,P,ay,f,e) , 


(3.10) 


with the four void density-profile parameters rs^Sc^a^f3^ the velocity-dispersion amplitude 
cr^, the logarithmic growth rate /, and the AP parameter s. For the width oj of the velocity 
dispersion in eq. (3.9) we use the best-fit values obtained from the data beforehand and do 
not vary it in the MCMC. We also do not distinguish between the velocity dispersion along 
or perpendicular to the line of sight, as the two are found to be very similar in fig. 4. Further, 
we assume a multivariate Gaussian likelihood function for the data given a model of the form 


= (27r) 


-^(Ivt-lvt)'c-i (ivt-evt) 


(3.11) 


where ^vt is the measured void-tracer cross-correlation function in redshift space, (^vt fhe 
corresponding model given by eq. (2.9), Ni^ its number of radial bins per dimension and C 
its covariance matrix. Note that the two-dimensional shape of ^vt h^^s been collapsed into 
the form of a vector in this notation, i.e. C is a matrix of dimensions x N^. We estimate 
the latter via Jackknife resampling of eq. (3.6), where the void stack is calculated Ny times 
after removing each individual void once. The resulting covariance matrix is then given by 


C 


A^v-1 

Ny 


E(^t-Ct) (c-Ct)', 


(3.12) 


where is the ith Jackknife sample and ^ Cvt mean of all Jackknife samples. 

Since C is only a statistical estimator of the true underlying covariance matrix C from the 
noisy data, it is subject to some degree of bias and uncertainty. The same is true for the 
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inverse of this matrix, which in some cases may not even exist. A convenient solution to this 
issue can be achieved via tapering of the covariance matrix, where its off-diagonal elements are 
increasingly diminished the more distant they are from the diagonal. An unbiased estimate 
of the inverse covariance matrix can be obtained by calculating [87, 88] 



where T is the so-called tapering matrix and o denotes an element-wise product. It can 
be shown that a specific type of functions yield a result that converges to the maximum- 
likelihood estimate of [87]. One such possible choice is the family of Wendland funetions 
[89], which in the case of d = 3 dimensions and /c = 1 reads 

'03,1 (^) = max (1 — r, 0)^ (1-h 4r) . (3.14) 

The elements of the tapering matrix can then be defined as Tij = 03p(|r^ — Yj \ /rt)^ where 
r = {r\\^r±y and rt denotes the tapering seale. Paz and Sanchez [ 88 ] have investigated 
the covariance tapering method on BAO measurements with galaxy auto-correlations, and 
found that it yields more accurate (and in particular lower) uncertainties in the model pa¬ 
rameters than the standard methods. In this case they report an optimal tapering scale of 
rt — 230h~^ Mpc, but point out that this value depends on the detailed structure of the 
covariance matrix. We have implemented the same value in our analysis and find it to yield 
reasonable results for void-galaxy cross-correlations, although it might not be the optimal 
choice. However, given this value to lie well within the linear regime of large-scale structure, 
where the two types of correlations are proportional to each other [67], this choice seems 
reasonable. The following uniform prior ranges for the values of our parameters are imposed, 

r, G [0.5,1.3]fv , 5c e [-1, -0.6] , a G [0, 3] , /3 G [5,13] , 

e [0,600]Jfi , / G [0,2] , £ G [0.7,1.3] . (3.15) 

The selection of these ranges is partly motivated by earlier numerical studies of the void- 
density profile from eq. (2.14) [62], and previous testing of our inference framework. 

Eventually, we explore the full parameter space with a Metropolis-Hastings algorithm 
provided in PyMC, and thereby sample the posterior probability distribution V{0\i^^t) of our 
model given the data with O(IO^) samples for each void stack separately. We also determine 
the maximum-likelihood model for each stack, which is represented by the white contours in 
figs. 1 , 2 and 3. As apparent from comparing the contour levels in these figures, the Gaussian 
streaming model yields an excellent fit to the data from the smallest to the largest voids in 
both dark matter and mock galaxies. In the following subsections, we will investigate the 
corresponding best-fit parameters and their uncertainties in more detail with the help of their 
full posterior distribution that resulted from our MCMC runs. 

3.5.1 Dark matter 

Let us begin by investigating the full posterior parameter distribution for one of our dark 
matter void stacks of radius range r^^ — (13.1 — 14.9)/i“^ Mpc shown in fig. 5. First of all, the 
void density-profile parameters (r^, 5c^ a, (5) are very consistent with what has been found in 
a dedicated one-dimensional analysis in ref. [62]. The constraints on these parameters are 
rather tight, showing no severe degeneracies amongst each other. Similarly, the posterior for 
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the cosmological parameter set (cr^, /, e) exhibits fairly symmetric and tight contours without 
strong degeneracies between any parameter pair. In particular, the preferred values for ^ 
3001^ are roughly consistent with what we measured directly from the velocities of the dark 
matter particles in our simulation as shown in fig. 4. Nevertheless, the logarithmic growth 
rate is not consistent with the expected value from general relativity, / = (z = 0) ~ 0.53, 
instead its posterior peaks at a value of zero. As already observed in fig. 1, the FoG effect 
dominates for void stacks in this radius range, leaving the Kaiser effect sub-dominant and 
hard to detect. On the other hand, the AP parameter more closely agrees with its fiducial 
value of e = 1 with a slight bias of about 6%. Hence, the 15% systematic error in the latter, 
as originally reported by Lavaux and Wandelt [14], can be reduced significantly when peculiar 
velocities are taken into account. Given the small relative Icr-uncertainty of only 0.5% in 
this measurement, this serves as a useful consistency check that our model is able to treat all 
anisotropic distortions accurately enough to allow conducting an AP test on the few-percent 
level. However, we realize that the accuracy provided by linear theory is not sufficient to 
yield unbiased parameter constraints of the model. 

Finally, it is interesting to note that the covariance between the two parameter sets 
discussed above is quite limited as well. This is expected, as the void density-profile param¬ 
eters in eq. (2.14) only influence the isotropic part of the redshift-space correlation function, 
whereas the cosmological parameters modify its anisotropic behavior. In principle, this al¬ 
lows to separate the inference process into two independent steps that can be iterated one 
after each other to find the optimal complete parameter set. 

The covariance matrix of the data is shown in the upper panel of fig. 8, normalized by 
its diagonal, which yields the correlation matrix Cij j ^CaCjy Because each two-dimensional 
void stack is represented by x ]S!\) bins, its covariance matrix exhibits x pixels, 
which can be thought of as a A'^ x N\) matrix of A"^ sub-matrices with the same dimensions. 
Since the covariance matrix is symmetric by definition, in the upper triangular part we plot 
the Jackknife estimate from eq. (3.12), while in the lower triangle we show its tapered version 
CijTij j^CiiCjj. It can be seen that the tapering has almost no effect on these scales, only for 
the very largest voids in our sample it significantly downsizes the far off-diagonal elements 
of the covariance matrix. The lower panel of fig. 8 depicts the inverse covariance matrix, 
again normalized by its diagonal. While the lower triangular part was obtained according 
to eq. (3.13), the upper part shows the direct inversion of eq. (3.12). Again, for the scales 
shown here and the high sampling density, the differences are marginal. 

3.5.2 Mock galaxies 

The posterior parameter distribution obtained from our dense mock-galaxy void stack with 
roughly the same range of effective radii — (13.0 — 14.8)/i“^ Mpc as for the dark matter 
voids above is presented in fig. 6. First of all we notice a very similar precision on all 
parameter constraints. Moreover, there is a slight decrease in r^, Jc, Oi and /3, resulting in 
a steeper void density-profile with a more pronounced ridge feature. This confirms earlier 
studies [84, 85] comparing void profiles in different tracer distributions. It can be explained 
by the biased clustering statistics of tracers of the mass distribution, resulting in an increased 
contrast in the density fluctuation of galaxies around voids. 

By examining the difference between figs. 5 and 6 one can also perceive an improve¬ 
ment in the inferred value of s, which is closer to its fiducial value in the dense mock-galaxy 
sample. While the constraining power on / is degraded, its maximum-likelihood value signif¬ 
icantly increases to a non-zero value. In fact, in the case of mock galaxies we are effectively 
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constraining the parameter ratio //6, which is often denoted j3 in the literature. This can 
directly be seen when plugging eq. (2.16) into eq. (2.12). The linear bias parameter of our 
dense mock-galaxy sample amounts to 6 0.8, which explains the relatively high value of the 

effective growth rate f /b ^ 0.67 here, a value which is perfectly consistent with the inferred 
posterior distribution within la. 

In contrast, the value of 6 1.8 for our sparse mock-galaxy catalog at higher redshift 

leads to a smaller effective growth rate of f/b = = 0.5)/b ^ 0.58, as can be seen 

in the posterior for the void radius range of = (48.1 — 53.1)h~^ Mpc in fig. 7. Although 
the high bias value leads to a steeper void-density profile, the lower galaxy number density 
and hence void abundance of this catalog results in wider parameter contours. Comparing 
the covariance matrices and their inverse in figs. 9 and 10, one can observe the stronger 
infiuence of the tapering procedure on the data from the sparse mock sample. Nevertheless, 
the posteriors for the logarithmic growth rate and the AP parameter nicely agree with their 
expected values within the 95.5% contour region. 

4 Discussion 

4.1 Joint constraints 

So far we have only considered the parameter inference from a single void stack within 
a narrow effective radius range. What remains to be investigated is how the parameter 
constraints among stacks from different void radii compare and whether it is possible to 
combine them in order to achieve a higher precision. 

Figures 11, 12 and 13 show the posterior distributions of / (respectively f/b) and 6: 
for each of our void stacks in dark matter, dense and sparse mock galaxies, respectively. 
Since these are the two main parameters of interest for the purpose of this paper, we do not 
investigate the remaining ones in detail anymore. As apparent from the different panels in 
these figures, the constraints on / and remain quite stable over the range of void sizes. 
While the Kaiser effect in dark matter voids can be only barely detected at the largest 
void radii, the growth rate inferred from voids in mock galaxies is quite consistent with its 
expected values from theory. For the sparse void sample, these values mostly fall within the 
95.5%-confidence region of the posterior distribution. The smallest voids are more prone to 
nonlinear effects, as they exhibit the highest surrounding ridges [62] that cause more extreme 
peculiar velocities (FoG) [82] . This can cause some discrepancy when using linear theory to 
model the void stacks in high-density tracers, as can be seen in fig. 12. We observe a slight 
bias of about 5% overestimation of the AP parameter and a varying bias underestimating the 
growth rate, depending on void radius. The best agreement is achieved for intermediate-size 
voids, which also provide the tightest constraints. It is a reassuring result in view of the fact 
that the preferred range of effective void radii for this type of analysis lies within the sample 
of voids we can identify. 

In order to combine the constraints from individual void stacks of a given tracer sample, 
we consider the joint likelihood as a product of the likelihoods of the individual stacks i, 

Aot=n'^(sv) • (4.1) 

i 

Note that we do not distinguish the model parameters 9 = f3,ay^ f^e) between 

different stacks. As argued above, the parameters of interest / and do not depend on void 
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size, so we can infer them from the joint dataset. This is of course not the case for the 
density-profile parameters (5)^ which all depend on void radius. However, we can 

simply marginalize over those, as we are not interested in their posterior distribution here. 

Figure 14 shows the joint posterior distribution for / and 6: from all void stacks obtained 
in our dark matter sample, as well as in our dense and sparse mock-galaxy catalogs. When 
compared to the constraints from the individual void stacks shown in figs. 11, 12 and 13, we 
observe a considerable gain in precision on the parameters of interest. In each of the mock- 
galaxy samples the la marginal error on f /b and 6: is reduced by roughly a factor of three. It 
amounts to a relative Icr-uncertainty of ^ 2% (20%) on f /b and ^ 0.2% (2%) on 6: (relative 
to their fiducial values) from the joint analysis of voids identified in the dense (sparse) mock- 
galaxy catalog. Thus, making use of the entirety of voids in a given catalog may significantly 
increase the information content that can be inferred via redshift-space distortions. However, 
using voids of all sizes in a single stack is suboptimal, as the redshift-space distortion signal 
from different physical scales is blurred out that way. 

4.2 Model details 

In order to examine the sensitivity to our model assumptions we have conducted a number 
of tests. First of all, we tried replacing the Gaussian form of the pairwise velocity probability 
distribution function in eq. (2.10) with an exponential distribution, as suggested by some 
authors (e.g. [43, 90-92]). Besides longer computing times in the MCMC sampling process, 
we found no noticeable differences in the final parameter constraints and concluded that 
prior knowledge on the exact shape of this function is of minor importance for the presented 
analysis. 

Similarly, we tested how sensitive our results are affected by the velocity dispersion 
model adopted in eq. (3.9), which depends on distance r from the void center, but ignores 
the direction dependence from eq. (2.17). By adding an additional parameter, we reran our 
MCMC analysis to constrain the two velocity dispersion amplitudes a\\ and a^ separately, 
and fixed the width uj of the Lorentzian in eq. (3.9) for each of its measured values along 
and perpendicular to the line of sight. In addition, we repeated the entire inference process 
with simply a constant cr^, neither dependent on scale, nor on direction. Apart from yielding 
slightly different best-fit values for respectively a\\ and cr^, these modification did not 
alter our parameter constraints in any notable fashion. This argues for the fact that the 
anisotropic shape of the void-tracer cross-correlation function is mainly determined through 
the mean of the pairwise velocity distribution function in eq. (2.10), and not its dispersion. 

Also note that the parameter space used in the presented inference process can be 
reduced. In ref. [62] it was found that the void density profile slopes a and [3 are strongly 
correlated with the scale radius for voids of different size. A simple parametrization 
from a fit to simulations was proposed to reduce the effective number of parameters by two. 
We implemented this approach into our MCMC analysis, but found almost no difference 
in the final constraints on neither /, nor 6:. This result is entirely expected, given the fact 
that we found no strong degeneracies between the void-density profile parameters and the 
cosmological parameters in sec. 3.5.1. 

5 Conclusion 

In this paper we presented a first attempt to analyze and model the effects of redshift-space 
distortions around cosmic voids with a focus on cosmological parameter inference. To this 
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end we generated catalogs of dark matter particles and realistic mock galaxies and used them 
to identify and stack voids in redshift space. In order to interpret the simulated mock data, 
we adopt the Gaussian streaming model for voids. We find that it successfully describes 
the void-tracer cross-correlation function in redshift space, given a realistic fitting function 
for the nonlinear real-space density profile of voids found in ref. [62], and its linear-theory 
relation to the corresponding velocity profile. All of these ingredients are then implemented 
into a Bayesian MCMC framework to infer the posterior parameter distribution of the model. 
The main conclusions drawn from this analysis are: 

• The values for the void density-profile parameters P) inferred from the two- 

dimensional void-tracer cross-correlation function in redshift space Cvt(f^||, are con¬ 
sistent with the values found in a dedicated one-dimensional analysis of the real-space 
void-density profile [62]. 

• Void density-profile parameters (r^, P) and cosmological parameters (cr^, /, e) show 
no strong degeneracies, because the former set describes the isotropic, the latter set the 
anisotropic part of <^vt(f ^||5 

• The growth rate / (respectively f /h) and the AP parameter 6: oc 1/DaH does not 
depend on void effective radius ry, in contrast to the void density-profile parameters. 
This allows to place joint constraints on these parameters from the entire available 
range of void sizes, yielding improvements in precision by factors of a few. 

• The relative uncertainties on f /b and 6: achievable in a survey volume of V = l/i“^Gpc^ 
with the presented mock galaxies range between /{f/b) ^ 0.02 — 0.2 and Ge/e = 
(^Dah/DaH 0.002 — 0.02. Increasing the survey volume, or adding independent in¬ 
formation from the CMB and the BAG, for example, will further reduce these numbers. 

• The 15% systematic error in the conventional AP analysis [14, 15] can be reduced to 
about 5% when peculiar velocities are taken into account. This type of measurement 
potentially contains a tremendous amount of information, consistent with earlier esti¬ 
mates [14]. However, the level of precision achievable, in particular for dense catalogs, 
falls below the model accuracy provided by linear theory. This argues for an extension 
of the current framework to the nonlinear regime with a more sophisticated model. 

It is important to note that the range of scales considered in this analysis is below the 
BAO scale, in a regime where theoretical models of the galaxy two-point correlation function 
fail (below 30h~^ Mpc to 50/i“^ Mpc, depending on redshift and galaxy bias, see e.g. [48- 
53]). We find that void-galaxy cross-correlations have the potential to unlock this scale 
range for precision cosmology. The principal advantages of void-galaxy cross-correlations 
are twofold: Firstly, even linear theory is able to describe the relation between density 
and velocity fluctuations around voids down to scales of ^ 10h~^ Mpc or lower [62], and 
secondly the contributions from redshift-space distortions only enter this statistic linearly, 
not squared, as the motion of void centers can be neglected [63]. In addition to that, the 
void auto-correlation function can be used to independently perform an AP test and further 
improve the constraints on 6: [18]. These advantages are, however, somewhat mitigated by 
the relatively low number statistics of voids compared to galaxies. 

Nevertheless, the analysis of void redshift-space distortions has the potential to con¬ 
tribute independent cosmological information from galaxy surveys that has so far remained 
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untapped. This method not only permits to improve our knowledge about the standard 
model of cosmology and gravity, but promises to test its alternatives. In particular, theo¬ 
ries of modified gravity predict deviations from general relativity to be most pronounced in 
unscreened low-density environments [93], making voids a smoking gun for the detection of 
fifth forces [94-101]. Redshift-space distortions of tracers around voids not only probe these 
low densities, they also contain information about the tracer velocities in that regime, which 
are directly related to the gravitational potential. Existing public void catalogs [102-106] 
already allow testing these scenarios, and future data from ongoing and upcoming galaxy 
surveys will extend the scope of this method. Combined with measurements of weak leasing 
to probe the underlying mass distribution inside and around voids [107, 108], such studies 
promise to shed more light on our concepts of dark energy, dark matter, and the laws of 
gravity. We will make our analysis tools publicly available by incorporating them into the 
VIDE pipeline on this website: https://bitbucket.org/cosmicvoids/vide_public/. 
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Figure 1. Dark matter void stacks Pvm (^||5 ^_L)/Pm — 1 in redshift space at 2 : = 0. Black solid/dashed 
lines show positive/negative contours of the data, white lines show the best-fit model. 
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Figure 2. Dense mock-galaxy void stacks Pvt{r\\,r±)/pt — 1 in redshift space at z = 0. Black 
solid/dashed lines show positive/negative contours of the data, white lines show the best-fit model. 
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Figure 3. Sparse mock-galaxy void stacks pwt{r\\^r±_)/pt — 1 in redshift space at 2 ; = 0.5. Black 
solid/dashed lines show positive/negative contours of the data, white lines show the best-fit model. 
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Figure 4. Velocity dispersion along (left) and perpendicular (right) to the line of sight around voids 
identified in dark matter particles at 2 : = 0 (top), dense mock galaxies at 2 : = 0 (middle), and sparse 
mock galaxies at 2 ; = 0.5 (bottom). Shaded regions depict the void-to-void standard deviation within 
each bin (scaled down by 20 for visibility), error bars show standard errors on the mean, and solid 
lines show the best fits from eq. (3.9). The insets list the mean effective void radius of each stack. 
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Figure 5. Posterior distribution of the model parameters for a dark matter void stack with effective 
radius range = (13.1 — 14.9)Mpc at 2 ; = 0. Solid, dashed and dotted contours show the 68.3%, 
95.5% and 99.7% confidence levels, respectively. Parameter values at the maximum of the marginal 
posterior distribution and its standard deviation are given at the top of each column. 
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Figure 6. Posterior distribution of the model parameters for a dense mock-galaxy void stack with 
effective radius range = (13.0 — 14.8)h“^ Mpc at z = 0. Solid, dashed and dotted contours show 
the 68.3%, 95.5% and 99.7% confidence levels, respectively. Parameter values at the maximum of the 
marginal posterior distribution and its standard deviation are given at the top of each column. 
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Figure 7. Posterior distribution of the model parameters for a sparse mock-galaxy void stack with 
effective radius range = (48.1 — 53.1)h“^ Mpc ai z = 0.5. Solid, dashed and dotted contours show 
the 68.3%, 95.5% and 99.7% confidence levels, respectively. Parameter values at the maximum of the 
marginal posterior distribution and its standard deviation are given at the top of each column. 


- 31 - 







































= (13.1—14.9)/i ^Mpc 


a * * « 

^ 4M4 

ii 


« # # 4 €^€jf€€€ 
# tf t 

j # 

jj«f— 
#»* 


W VV VVV-.r^^V'tVV 
fW^r'r V-'r-'. VVV^r-'r'^. W^rV 


V-'r V V V ■‘r V'',. V W VV' 
, ■V'‘r-V-V'VVV'.'.'tV< 

V W • 

VVV^r^VV^v-*-^-- 

W ' r-'.r ■<, ^ VV'-r 

“'X'V'. •■ 

0_ i. ■'. .•- -L ■- V 'r M- V*' 


- ■ vv'r'‘r<yvv-- ■' ■--. 


•'. '. • VV-'r > 

. i i J 4^ J i 



Q 


m 

I 

m 


k 


mv ///// 


v v v v -r v 

r»V''' ■*-'- ‘- *-*- 

jfc" 

■ * ■ ' 


0.0 


0 5 10 15 20 25 

rjL [/i“^Mpc] 


= (13.1—14.9)/i ^Mpc 




iO.5 

0.3 

0.2 

^0.0 

- 0.2 


10 

r^ [/i“^Mpc] 


0.9 

0.8 

0.6 

0.5 

0.3 

0.2 


Figure 8. Covariance matrix (top) and its inverse (bottom) of a dark matter void stack with radius 
range r^ — (13.1 — 14.9)h“^ Mpc at z = 0. Each matrix is normalized by its diagonal, the lower 
triangular parts have been obtained using the tapering technique as described in the text. 
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Figure 9. Covariance matrix (top) and its inverse (bottom) of a dense mock-galaxy void stack with 
radius range = (13.0 — 14.8)h“^ Mpc at 2 : = 0. Each matrix is normalized by its diagonal, the 
lower triangular parts have been obtained using the tapering technique as described in the text. 
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Figure 10. Covariance matrix (top) and its inverse (bottom) of a sparse mock-galaxy void stack 
with radius range = (48.1 — b3.1)h~^ Mpc at ^ = 0.5. Each matrix is normalized by its diagonal, 
the lower triangular parts have been obtained using the tapering technique as described in the text. 
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Figure 11. Constraints on growth rate / and AP parameter 5 for each void stack in the dark matter. 
White stars show fiducial values of / = = 0) = 0.53 and 5 = 1. Confidence regions as shown 

only include statistical uncertainties and marginalize over the shape of the void density profile, the 
residual offset is due to modeling error (as discussed in the text). 
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Figure 12. Constraints on growth rate f/b and AP parameter 5 for each void stack in the dense 
mock-galaxy sample. White stars show fiducial values of f/b = = 0)/0.8 = 0.67 and 5 = 1. 

Confidence regions as shown only include statistical uncertainties and marginalize over the shape of 
the void density profile, the residual offset is due to modeling error (as discussed in the text). 
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Figure 13. Constraints on growth rate f/b and AP parameter 5 for each void stack in the sparse 
mock-galaxy sample. White stars show fiducial values of f/b = = 0.5)/1.8 = 0.58 and 5 = 1. 

Confidence regions as shown only include statistical uncertainties and marginalize over the shape of 
the void density profile, the residual offset is due to modeling error (as discussed in the text). 
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Figure 14. Joint constraints on growth rate / (respectively f/b) and AP parameter e for all void 
stacks from dark matter at z = 0 (top), dense mock galaxies at 2 ; = 0 (bottom left), and sparse 
mock galaxies at ^ = 0.5 (bottom right). Parameter values at the maximum of the marginal posterior 
distribution and its standard deviation are given at the top of each column. 
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